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We continue our earlier study of the phase structure of a SU{2) gauge theory whose action 
contains additional chirally invariant four fermion interactions. Our lattice theory uses a reduced 
staggered fermion formalism to generate two Dirac flavors in the continuum limit. In the current 
study we have tried to reduce lattice spacing and taste breaking effects by using an improved fermion 
action incorporating stout smeared links. As in our earlier study we observe two regimes; for weak 
gauge coupling the chiral condensate behaves as an order parameter differentiating a phase at small 
four fermi coupling where the condensate vanishes from a phase at strong four fermi coupling in 
which chiral symmetry is spontaneously broken. This picture changes qualitatively when the gauge 
coupling is strong enough to cause confinement; in this case we observe a first order phase transition 
for some critical value of the four fermi coupling associated with a strong enhancement of the chiral 
condensate. We observe that this critical four fermi coupling varies monotonically with bare gauge 
coupling - decreasing, as expected, as the gauge coupling is increased. We have checked that these 
results remain stable under differing levels of smearing. They appear to rule out the appearance of 
new fixed points associated with chirally invariant four fermion interactions in confining non abelian 
gauge theories. 

PACS numbers: 05.50.+q, 64.60.A-, 05.70.Fh 



I. INTRODUCTION 



The primary purpose of the Large Hadron Collider 
(LHC) is to probe the nature of electroweak symmetry 
breaking in the Standard Model (SM). In the SM this 
is accomplished using a single scalar field - the Higgs. 
Unfortunately scalar field theories are known to suffer 
from fine tuning and triviality problems and in conse- 
quence other extensions of the SM based on supersym- 
metry or new strong dynamics have been put forth to pro- 
vide more natural explanations of electroweak symmetry 
breaking. One of the simplest realizations of dynamical 
electroweak symmetry breaking are technicolor models 
in which the SM is augmented with new fermions (tech- 
niquarks) which interact via a new strong technicolor 
gauge force. Strong interaction effects can then drive 
the formation of a techniquark chiral condensate which 
will in turn precipitate electroweak symmetry breaking 
if the electroweak gauge group appears as a subgroup of 
the chiral group of the techniquarks [TJ [2] . While theo- 
ries with Q CD-like dynamics are experimentally excluded 
they are far less stringent constraints on these models if 
the dynamics deviates dramatically from QCD. This is 
particularly true in light of the recent discovery of a light, 
Standard Model-like Higgs at the LHC. One way to pro- 
duce such a light state in a strongly coupled BSM theory 
is to consider models lying close to the lower boundary 
of the conformal window. In this case it is conjectured 
that one might be able to identify the experimentally ob- 
served Higgs with a pseudo Goldstone boson - the dilaton 
- arising from spontaneous breaking of approximate con- 
formal symmetry. Such theories are often termed walk- 
ing gauge theories since the gauge coupling should evolve 
slowly with energy scale [TJ [3] and possess other poten- 
tially useful features such as reduced contributions to the 



S parameter and flavor changing neutral currents. 

The search for such walking gauge theories has led to 
numerous recent lattice studies - see the conference re- 
views [3H7]. However, the simplest technicolor models, 
while giving rise to masses for the gauge bosons, are in- 
capable of generating masses for SM fermions. To achieve 
this more elaborate schemes such as extended technicolor 
(ETC) [8HTT] and top-condensation [T2HT5] have been 
invoked. At low energies these latter theories share a 
common framework; the actions contain additional four 
fermion interactions coupling both SM fermions and tech- 
niquarks. Indeed most theories of BSM physics can be 
expected to yield effective actions at low energy which 
contain additional four fermion terms. 

In light of this it is important to study the structure 
of non abelian theories in the presence of additional chi- 
rally invariant four fermion interactions. This was the 
motivation behind our earlier work [16] . where we ex- 
plored the phase structure of a NJL model gauged under 
the SU (2) group. This model has also been examined us- 
ing truncated Schwinger-Dyson calculations in [17] where 
strong four fermi interactions were claimed to lead to 
new universality classes of model with potentially inter- 
esting consequences for model building. Our work can 
also be seen as a non-perturbative check on these trun- 
cated Schwinger-Dyson results. 

In our previous work we observed that the chiral con- 
densate of the strongly coupled gauge theory underwent a 
strong enhancement for large values of the four fermi cou- 
pling. The transition appeared very abrupt and therefore 
likely first order but we were unable to find definitive evi- 
dence of this in the earlier work. In addition we observed 
that the critical four fermi coupling did not vary mono- 
tonically with the gauge coupling which led us to suspect 
that the results of those simulations might be plagued 
with large lattice spacing effects. To counter this we have 
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revisited the model using an improved fermion action in 
which the usual gauge links are replaced by so called stout 
smeared gauge links. The smeared gauge fields are ob- 
tained by a gauge invariant averaging procedure carried 
out over a local region in space. The essential idea is 
that this averaging procedure will tend to remove large 
UV fluctuations from the lattice gauge field while pre- 
serving the IR physics. Smeared actions been shown to 
be very effective in lattice QCD [T9H2T] . The specific 
scheme we employ was developed in [22] for the case of 
5/7(3). 

In the following section, we briefly describe our lat- 
tice theory and its connection to the gauged NJL model. 
Next we describe the method of stout smearing in the 
context of dynamical fermion simulations of an SU(2) 
lattice gauge theory. Finally, we move on to our numeri- 
cal results, where we discuss the phase structure resulting 
from varying both gauge and four fermi couplings. We 
show that suppressing lattice spacing effects via smear- 
ing removes some of the technical issues associated with 
our earlier work but does not change the final conclusion: 
that the non abelian theories do not develop new fixed 
points associated with strong four fermi interactions. In- 
stead, the lattice theory undergoes a first order phase 
transition for large four fermi coupling. This transition 
is associated with a rapid enhancement of the chiral con- 
densate and the generation of large fermion masses. 



II. FOUR FERMION INTERACTIONS 
AND OFF THE LATTICE 



ON 



We will consider a model which consists of two flavors 
of massless Dirac fermions transforming in the funda- 
mental representation of a SU (2) gauge group and cou- 
pled through an SU (2) l x SU (2)r chirally invariant four 
fermi interaction of NJL type. To avoid potential sign 
problems we have utilized two copies of this basic 2 fla- 
vor doublet yielding a four flavor theory. The action for 
a single doublet takes the form 



S = I d 4 x - M - 



1 

2g 



2 Tr[F^FH, 



(1) 



where G is the four-fermi coupling, g the usual gauge 
coupling and r a ,a = 1...3 are the generators of the 
SU{2) flavor group. 

To implement the theory in Eq. [I] on the lattice, it is 
convenient to re-parametrize the four fermi term in the 
continuum action via the use of scalar auxiliary fields. 
Specifically the fermion interaction term is replaced by 
the Yukawa terms 



G 



(Wfa + &lsT a ilxt> a )+- {<t>\ + <t>; 



(2) 



We discretize the above action using a (reduced) stag- 
gered fermion formalism. This approach allows build a 
lattice theory which admits just two continuum flavors 
of fermion in the continuum limit without encounter- 
ing rooting issues. Furthermore it allows us to write 
down lattice four fermi terms which are invariant un- 
der a discrete subgroup of the continuum chiral symme- 
tries [23 , 24 . The presence of four fermi interactions has 
an additional attractive feature - it allows us to study 
the lattice theory with exactly zero fermion mass where 
these discrete symmetries are exact [25] . 

In order to do this, we first rewrite the fermionic sector 
of the continuum theory in terms of a new set of matrix 
valued fields, ^ and ^ which naturally arise after twisting 
the original Lorentz symmetry with the flavor symmetry. 
We then expand these matrices on a basis correspond- 
ing to products of gamma matrices and associate these 
products with staggered fields X- 

y(x) = lY^ x+b x(x + b), (3) 
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(4) 



where j x + b = Yl i=1 ^ i+hi and the sums correspond to 
the vertices in an elementary hypercube associated with 
lattice site x as the components vary b{ = 0,1. The 
lattice theory thus far describes four flavors of continuum 
fermion. To reduce the flavor content further it is possible 
for massless fermions to restrict the single component 
fields x an d X t° even and odd sites, respectively. In this 
way we arrive at the reduced staggered fermion action: 



(5) 



where, 
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b) 



(7) 



and rjn(x) is the usual staggered quark phase given by 



(8) 



Notice that the action involves a single staggered field x 
defined over all lattice sites but coupled through a modi- 
fied gauge field U fI (x). The two staggered tastes become 
the two physical quark flavors in the continuum limit. 
This action is invariant under the discrete shift transfor- 
mations: 



X(x) £p(x)x(% + p), 
U^x)^U;(x + p), 



(9) 
(10) 

(11) 
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These shift symmetries correspond to a discrete subgroup 
of the continuum axial flavor transformations which act 
on the matrix field \I/ according to 

* (12) 

III. STOUT SMEARING 

To reduce the magnitude of lattice spacing effects we 
replace the links (x) used in constructing the staggered 
fermion action by stout smeared links U'^{x) which are 
defined as 



The parameter p can be tuned to optimize the smearing 
- in this work we have set p = 0.1. To leading order in p 
it is easily seen that this procedure generates a smeared 
link which is a weighted sum over the original link and its 
surrounding staples. In practice this procedure is iterated 
with the level n smeared fields being constructed 

from the fields ujf 1-1 " 1 at level n-1 via the obvious relation 

U^\x) = {e^^^jU^ix) (17) 

We have run our simulations using n = 1 and n = 2 lev- 
els of smearing. For details on how this smearing in im- 
plemented and the calculation of the smeared force term 
required by the RHMC simulation we refer the reader to 
the Appendix. 

IV. NUMERICAL RESULTS 

Upon integration over the basic fermion doublet we ob- 
tain a Pfafflan Pf(M(U)) depending in a complicated way 
on the original gauge field L^(xrj The required pseud- 
ofermion weight for 2 such doublets (Nf = 4 flavors) is 
then Pf(M) 2 . The pseudoreal character of SU(2) allows 
us to show that this Pfafflan is purely real and so the four 
flavor simulations suffer from no sign problem. In prac- 
tice we implement this weight in the path integral using 
a pseudofermion operator of the form (M^M)~^ and a 
RHMC algorithm. A standard Wilson gauge action is 
employed for the gauge fields. 



Note that the fermion operator appearing in eqn. |5] is antisym- 
metric and hence a Pfaffian rather than a determinant results 
after integrating over the single fermion field x 



where the traceless, antihermitian matrix is given by 



QM = \ pi], (x) - Q^x)] - ± Tr[Qt ( X ) - %(x)]. 

(14) 

with 



Stv(x) = Cp(x)Ut(x), (15) 



(16) 



I 

We have focused our study on a lattice of size 6 4 and 
a bare fermion mass m = 0. We have simulated a range 
of gauge couplings f3 = 4/g 2 = 1.8 — 2.2 with additional 
runs at f) = 10.0 to test the pure NJL model in which the 
gauge interactions are switched off. For each gauge cou- 
pling we have performed an extensive scan over a dozen 
or so values of the four fermi coupling G. The range of 

G = 0.1; Nf = 4 




FIG. 1. Polyakov line vs /3 at G = 0.1 for four flavours 

/3 that was employed in the simulations was determined 
by examining the average plaquette and Polyakov line 
as f3 was varied holding the four fermi coupling fixed at 
G = 0.1. In figure [l] we see a strong finite volume 
phase transition for f3 > 2.2 on a lattice with L — 6 
and no smearing. Furthermore, this transition occurs at 
smaller bare gauge coupling, j3 > 1.9, when two levels of 
stout smearing are employed. 

The primary observable used in this study is the chiral 



K( x ) = {z Q ' {x) }U^x). (13) and 
I 

C m = P U ^ x + v ) U »( x + V) + Ut{x - v) U^x - v) U v {x -v + p% 
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condensate which is computed from the gauge invariant 



6 4 ;/3 = 10.0 



1.5 " 




FIG. 2. |(xx)| vs G for p = 10.0, 6 4 lattice 



Because of the absence of spontaneous symmetry break- 
ing in finite volume we measure the absolute value of this 
operator. In a chirally broken phase we expect this to ap- 
proach a constant as the lattice volume is sent to infinity. 
Conversely if chiral symmetry is restored this observable 
will approach zero in the same limit. We have also used 
the auxiliary scalar field to monitor this observable. 

Having used these simple gauge observables to deter- 
mine an appropriate range for the gauge coupling f3 we 
now turn to the behavior of the system as the four fermi 
coupling is varied. As a benchmark we have examined 
the system for a small value of the gauge coupling at 
j3 = 10.0 The condensate is shown in figure [2] as the 
four fermi coupling varies. The plot reveals evidence for 
a phase transition separating a chirally symmetric phase 
at small G from a broken phase for G > G c . The rel- 
atively smooth behavior is consistent with earlier work 
using sixteen flavors of naive fermion reported in [26] 
which identified a line of second order phase transitions 
in this region of parameter space. It also agrees with the 
behavior seen in previous simulations of staggered quarks 

We now turn to our results at non-zero gauge coupling, 
corresponding to smaller values of /3 (stronger gauge cou- 



2 Note: gauge coupling, g 2 varies as i. 



one link operator 



( t( x ~ e M )x(^ - e M ))> e(x)^(x) (18) 
I 



6 4 ;f3 = 1.8 
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FIG. 3. (xx> vs G at = 1.8 for 6 4 lattice 



pling). Figure [3] shows results for the condensate in the 
confining theory at /3 = 1.8 for several levels of smear- 
ing. Notice that the chiral condensate is now non-zero 
even for vanishing four fermi coupling consistent with 
spontaneous chiral symmetry breaking in the pure gauge 
theory. However, it jumps abruptly to much larger values 
when the four fermi coupling exceeds some critical value. 
Notice that while the limiting values of the condensate 
in both small G and large G phases do not depend on 
how much smearing is done, the position of the transi- 
tion itself depends on the smearing level. This should 
not be particularly surprising; the magnitude of critical 
couplings associated with lattice phase transitions gener- 
ically depend on details at the scale of the lattice spacing. 
However, notice that the chiral condensate, being an IR 
quantity, does not show this sensitivity. 

Figure [3] appears to show that the jump in the chiral 
condensate is markedly discontinuous in character - remi- 
niscent of a first order phase transition. Further evidence 
for this first-order behavior is seen in figures |4|5[ |6j and 
[7] which show histograms of the condensate (as measured 
by the scalar field average a 4) for several four fermi cou- 
plings G around the phase transition. The histograms 
tell a clear story; for G both below and above the transi- 
tion a single peak is seen. In contrast for values of G close 
to the jump in the condensate the double peak structure 
indicates that condensate has tunneled back and forth 
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6 4 ;/? = 1.8; G = 0.72 



6 4 ; p = 1.8; 67 = 0.80 





(2 — level smearing) 



-0.03 0.13 

(2 — level smearing) 



FIG. 4. Histogram of a 4 for G=0.72. 



FIG. 6. Histogram of a 4 for G=0.80. 



6 4 ;/3 = 1.8; G = 0.76 



1.8; G= 1.44 





cr 4 

(2 - level s 



-0.55 -0.5 

(2 — level smearing) 



FIG. 5. Histogram of a 4 for G=0.76. 



FIG. 7. Histogram of a 4 for G=1.44. 



between two different states over the course of the run. 
Such a tunneling is clear indication of the presence of a 
first-order transition in the lattice system. 

We have also examined the dependence of these fea- 
tures on the magnitude of the gauge coupling f3. Fig- 
ure [8j shows a plot of the condensate vs G for Nf = 4 for 
a fixed lattice size L = 6 and using two levels of smear- 
ing. Two regimes are clearly visible; for j3 < 1.9 or so 
the chiral condensate is non zero for small G correspond- 
ing to a strongly coupled, confining and chirally broken 
gauge theory even in the absence of four fermi terms. 
Conversely for f3 > 2.0 the chiral condensate goes to zero 



for vanishing four fermi coupling and the system is de- 
confined in the absence of the four fermi coupling. This is 
not a surprise; it corresponds to the usual deconfinement 
one sees when the lattice size becomes smaller than the 
confinement scale. 

Notice that the onset of the crossover occurs for smaller 
values of the critical four- fermi coupling, G cr , as we in- 
crease the strength of the gauge coupling (corresponding 
to decreasing f3) . This monotonic dependence of the crit- 
ical four fermi coupling contrasts with the situation seen 
earlier in our study with an unsmeared action [16]. In 
the strong coupling regime, since the chiral symmetry is 
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6 4 ; 2-level smearing 




FIG. 8. (xx) vs G for varying /3. 2-level smearing has been 
implemented here. 



already broken, the cross over occurs quickly at a small 
(weak) four-fermi coupling. As one approaches the weak 
coupling limit, the model is less influenced by the gauge 
interactions and the chiral symmetry remains intact until 
the four-fermi coupling becomes strong enough to break 
the symmetry. 



V. SUMMARY 

In this paper we have conducted numerical simulations 
of a four flavor non abelian gauge theory in the pres- 
ence of additional chirally invariant strong four fermion 
interactions. The latter are expected to be present in 
generic effective theories of BSM physics including (walk- 
ing) technicolor, top quark condensation and little Higgs 
theories. We employ a reduced staggered fermion lattice 
action which preserves a discrete subgroup of the contin- 
uum chiral symmetry. In addition we have improved this 
action using stout smearing techniques. 

We have examined the model for a modest lattice size 
of 6 4 and over a range of gauge couplings /?, and four 
fermi interaction strength G. In the NJL limit f3 — >• oo 
we find evidence for a continuous phase transition cor- 
responding to the expected spontaneous breaking of chi- 
ral symmetry due to strong four fermi interactions. The 
location and character of this transition is, not surpris- 
ingly, insensitive to the smearing. However, for gauge 
couplings (3 that generate confinement and a non-zero 
chiral condensate even at G — 0, this transition appears 
much sharper and indeed we provide strong evidence that 
the transition is first order. Such a phase transition pre- 
cludes the existence of any new continuum limits in the 
model. The location of this transition occurs at a critical 
four-fermi coupling which depends on the level of smear- 



ing and decreases monotonically with increasing gauge 
coupling. 

Thus our results are consistent with the idea that the 
second order phase transition which exists in the pure 
NJL theory (G = oo) survives in the presence of weak 
gauge interactions. However our results indicate that this 
line of transitions becomes first order at a point at which 
the gauge coupling becomes strong enough to cause con- 
finement. In a finite volume this occurs at a finite /? but 
presumably this first order region extends all the way to 
/3 = = oo in infinite volumes. The gauged NJL model 
we study in this paper was examined using a truncated 
Schwinger-Dyson analysis in [17] in which it was claimed 
that the phase diagram of the gauged NJL model should 
include a new critical line along which critical exponents 
such as the mass anomalous dimension vary continuously. 
Our results are inconsistent with this scenario but are 
consistent with the approximate analysis performed in 
[29] which argues for a strong enhancement of the chi- 
ral condensate for large four fermi coupling. While our 
results seem to exclude new universality classes due to 
strong four fermi couplings in confining gauge theories 
it is still logically possible that this conclusion might be 
very different in a gauge theory which is conformal at 
G = 0. 
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Appendix A: Calculation of the stout smeared force 

in SU{2) 

Stout smearing is implemented in our code in two 
parts. The first involves replacing the gauge links in the 
fermionic action by the smeared links, ujP\x) as given 
in equation [IT] The second, and the non-trivial, part 
of implementing smearing lies in the calculation of the 
smeared force F^(x) needed in the auxiliary molecular 
dynamics evolution which forms the kernel of the dy- 
namical fermion algorithm. The force F fI (x) is defined 
as: 

f ^ = u ^ 9 -wm (A1) 

The complication that arises at this point is that we have 
no explicit representation of the smeared action in terms 
of the original gauge links - indeed such an expression 
would be very complicated given that the smeared links 
are determined iteratively. Instead we follows the ap- 
proach of [22 and derive an iterative expression for the 
smeared force by considering the rate of change of the 
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action over a fictitious (simulation) time t: 



The matrix exponentials can be handled using the SU(2) 
identities: 



sing 



(cosq)I 2 x2 H Q 



(A4) 



at 



S(U (k) ) = 2 Re^ Tr 



(A2) 



where q the eigenvalues of the traceless antihermitian 
matrix Q are given by 



1 



Tr(Q 2 ) 



(A5) 



Here, $\x) 



as(u {k) ) 



w ~ In ° rder t0 Simplify notation ' Using this the derivatives of e« can be expressed 

we denote the k level fields with a prime and the {k — 1) 
level fields without a prime. Using equation ( [l3| ), we 
have, 



d 



at 

with 



dt 



B 



as 

sinq / dQn(x) 



q 



dt 



a_ 
at 



_ sin_9 _ (cos(/ - sin (//</) 



(A3) 



2g 2<? 2 
Using these results eqn. can be expressed as 



(A6) 
(A7) 



£ (/iWe^^l^w) - Tr (a 



at 



(A8) 



where 

A M (ar) = - [rt ( X ) - T^x)] 
and 



N 



Tr [rt ( X ) - rv(x)] , 

(A9) 



r M (ar) = Tr [f^B^U^x)} Q^x) 



(A10) 

and ^(x) is given by equation (15). Using the cyclic 
properties of the trace we ultimately obtain, 



fM = fl{x){cosqI 2X2 + S ^Q) - Cl(x)T^x) 



p^2e^(x) 

jJL,X 



where, 



+Ul(x + /i - v)Ul(x - v)K{x - v)U^(x - v) 
-Ul(x -fJ> + v)Ul(x - v)K{x - v)XJ v (x - v) 
-k v [x + \i)U v (x + fi)Ul(x + ^)^(x) 
+ + v)A^(x + ^)^(x) 

+ /x - ^)A„(x + /i - v)Ul(x - v)U v {x - v). 

(A12) 



All 



gives f^x) = in terms of f' (x) 



Eqn. 

ffl k \x) at level Given an initial "naive" force at level 

n = dS<yU (t) } it may be iterated to yield (x) and 

hence the final smeared force needed by the molecular 
(AH) 

dynamics equations 



F^x) = U^x)fj?\x), 



(A13) 
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